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The behavior of the quintessence field is studied during inflation. In order to have a satisfactory 
model of dark energy, the quintessence field value today should be as insensible to the initial con- 
ditions as possible. Usually, only the dependence on the initial conditions specified at the end of 
inflation or, equivalently, at the beginning of the radiation dominated era, is considered. Provided 
the quintessence field is initially within a large but, crucially, finite interval, its present value becomes 
independent of the initial value it started from. The question as to whether inflation naturally drives 
the quintessence field to the above-mentioned interval is addressed. Since the quantum effects turn 
out to be important, the formalism of stochastic inflation is used in order to calculate the evolution 
of the quintessence field. Moreover, the quantum effects originating from the inflaton field are also 
taken into account and are proved to be sub-dominant in most cases. Finally, the requirement that 
the quintessence field is on tracks today is shown to imply quite tight constraints on the initial 
values of the quintessence and inflaton fields at the beginning of inflation. In particular, the initial 
value of the inflaton field cannot be too large which indicates that the quintessential scenario seems 
to be compatible with inflation only if the total number of e-folds is quite small. 

PACS numbers: 98.80.Cq, 98.70.Vc 



I. INTRODUCTION 

The observations suggesting that our Universe is 
presently undergoing a phase of accelerated expansion 
have recently accumulated Q,H[13- If really confirmed, 
this discovery is certainly a breakthrough for cosmology 
but, at the same time, represents a big challenge since 
finding a convincing explanation for such a phenomenon 
is clearly a difficult task. 

From a theoretical point of view, the presence of a 
non-vanishing cosmological constant whose energy den- 
sity would be of the order of the critical energy density 
today seems to be the most natural solution. In addition, 
the currently available data on the equation of state are, 
so far, compatible with this assumption. But it is well- 
known that the theoretical preferred value of the cosmo- 
logical constant corresponds to an energy density much 
larger than the critical energy density and there exists, 
at the moment, no convincing arguments which would 
explain this difference 

This situation has led the physicists to seek for alter- 
natives. Among the solutions proposed, the quintessence 
scenario has recently attracted a lot of attention [g, LJ, B 
H El El El El El It consists in postulating that the 
acceleration of the expansion is caused by a scalar field, 
the quintessence field Q, evolving in a potential the typi- 
cal shape of which is given by W(Q) = M 4+a Q~ a , where 
M is an energy scale and a > a free parameter |6] . The 
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main advantage of this scenario is that the coincidence 
problem can be solved because the equations of motion 
possess a solution which is an attractor. Therefore, the 
present evolution of the quintessence field is independent 
from the choice of the initial conditions. Moreover, when 
the field is on tracks, Q is typically of the order of the 
Planck mass and hence, for not too small values of a, 
the scale M can be large. As a result, the fine-tuning is 
less severe than in other scenarios because it is possible 
to explain the presence of a very small scale (the vacuum 
energy density today) by means of a theory which, on 
the contrary, is characterized by a large scale M. This 
is due to the inverse power-law shape of the potential 
and is reminiscent of the "see-saw" mechanism in par- 
ticle physics. This has also the advantage that model 
building can be considered in the realm of high energy 

physics H El El El El Ip- 
so far most of the studies have been devoted to un- 
derstanding how the quintessence field evolves from the 
beginning of the radiation era until now. Another impor- 
tant (related) question, in view of its observational im- 
plications, has been to estimate the value of the equation 
of state today. In this paper, we address a new question, 
namely that of the behavior of Q during (chaotic) infla- 
tion assuming, for simplicity, that the quintessence field 
and the inflaton are not coupled. This is an important 
problem since it is crucial to check that Q is, at the end 
of inflation (or at the beginning of the radiation domi- 
nated epoch), in the range of values which are such that 
the field is on tracks today. 

However, the problem does not only boil down to solv- 
ing the Klein-Gordon e quat ion in an inflationary back- 
ground. Indeed, in Ref. 15], it has been suggested that 
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the quantum effects could play an important role. In this 
case, the techniques of stochastic inflation can be 

used to describe the evolution of the quin tessence field. 
This was done for the first time in Ref. . The method 
utilized in that article was to solve the Fokker-Planck 
equation in order to follow the evolution of the proba- 
bility distribution of the quintessence field, ft was then 
shown that, typically, the attractor is joined at relatively 
small redshifts. 

In the present paper, we consider the above-mentioned 
question again but from a different perspective. One of 
our main purposes is to calculate the probability distri- 
bution function of the quintessence field at the end of 
inflation (or at the beginning of the radiation dominated 
era) in order to estimate whether it is likely that the value 
of Q is such that the attractor is joined today. Moreover, 
requiring that the corresponding probability is significant 
can be used to constrain the space of the initial condi- 
tions, i.e. the initial value of the inflaton (or, equivalcntly, 
the total number of e-folds) and quintessence fields. In 
addition, we demonstrate that this also puts constraints 
on the power index a characterizing the shape of the 
quintessence potential, namely small values of a are dis- 
favored. 

Another goal of the present work is to include the in- 
flaton fluctuations, to study under which circumstances 
their effect can be important and, when it is the case 
(and when it is possible), to calculate the correspond- 
ing correction to the behavior of the quintessence field. 
Indeed, in Ref. ^^|, since the inflaton was treated as 
a classical field, the quantum effects were only sourced 
by the quintessence noise. However, the inflaton itself 
is also influenced by the quantum effects and, there- 
fore, a priori the inflaton noise also affects the evolution 
of the quintessence field. In fact, the variance of the 
quintessence field can be written as 

<r\t) = -^j\*[<p{T)]&T, (1) 

where H is the Hubble parameter. Roughly speaking, 
taking into account the inflaton noise amounts to put the 
coarse-grained inflaton in the above equation rather that 
its classical counterpart. If the expression calculated in 
this way differs significantly from the expression obtained 
by inserting the classical inflaton, then the inflaton noise 
plays indeed a non negligible role. 

Another difference from Ref. is that we directly 
solve the Langevin equation rather than the Fokker- 
Planck equation. Obviously, this is only a technical dif- 



ference since the two approaches are equivalent. In a first 
time, the Langevin equation is solved by means of a per- 
turbative expansion. In this regime, we show that the 
influence of the inflaton noise is always negligible. Then, 
in a second time, we try to solve the Langevin equation in 
the non-perturbative regime (for the quintessence field) 
by modeling the effect of the classical force with a wall. 

This article is organized as follows. In the next section, 
we quickly review the basic principles and equations of 
the stochastic approach. Then, in Sec. Ill, we present the 
perturbative method used to solve the Langevin equa- 
tion. We apply this method to inflation, compare the re- 
sults obtained with those already known in the literature 
and demonstrate that they are equivalent. In Sec. IV, we 
apply the perturbative approach to the Langevin equa- 
tion for the quintessence field and explicitly calculate the 
quintessential quantum effects. In this regime, we show 
that the contribution coming from the inflaton noise is 
negligible. Then, we present a model with a reflecting 
wall which allows us to explore a region where the pertur- 
bative approach breaks down. We study the constraints 
on the initial conditions of the inflaton and quintessence 
fields that exist in order for the coincidence problem to 
still be solved. We prove that these constraints are quite 
stringent. Finally, in Sec. V, we discuss the results ob- 
tained in this article and present our conclusions. 



II. BASIC EQUATIONS 

In the Friedman-Lemaitre-Robertson- Walker (FLRW) 
Universe, the metric of which can be written as ds 2 = 
— dt 2 + a 2 (t)5ijdx' t dx : > (we assume flat space-like sec- 
tions), the evolution of a scalar field </>(£, x) is described 
by the Klein-Gordon equation 

<p + 3H(f) - — I — — - = ) (2) 

a z dip 

where a dot denotes the derivation with respect to the 
cosmic time t. 

In the stochastic formalism ^3]> one is interested in 
the dynamics of a "coarse-grained" field (p(t,x). This 
coarse-grained field is defined to be the average of the 
ordinary field <j) over a physical volume whose size is 
somewhat larger than the Hubble radius H^ 1 = a/ a. 
Therefore, ip(t, x) basically contains the long-wavelength 
Fourier modes (with wavenumber k < aH) only. Techni- 
cally, one writes 



(2tt) 3 / 2 



dkW (k - oaH) [c kA i fc (i)e !k x + 4/4 (t)e" 



ik-x 



(3) 



where W(z) is the so-called window function. In the case of a white noise, the window function is the step function. 



3 



In a more realistic situation, the window function should 
be taken as a smoothed version of the step function . 
This corresponds to the case of a colored noise and the 
problem is generally technically more complicated in this 
situation. In this article, for simplicity, we restrict our- 
selves to the case of a white noise. It should also be 
noticed that, in Eq. @, the mode functions //fe(t) are, by 
definition, the free mode functions, i.e. obey the equation 
flk + 3i?/ifc + (k 2 /a 2 )fik — 0. Finally, a is a parameter 
smaller than 1 , introduced in order to allow some level of 
arbitrariness in the choice of the smoothing scale. 

The evolution of the coarse-grained field is still de- 
scribed by the Klein-Gordon equation but a suitable 
random noise field acting as a classical stochastic 
source term, should be added to the right hand side in 
order to mimic the quantum fluctuations. In the slow- 
roll approximation (p is negligible compared to 3H (p and, 
since we are dealing with super-Hubble scales, the gradi- 
ent term can also be dropped. The coarse-grained field is 
thus governed by a first order Langevin-like differential 
equation which can be put in the form 



d^ 1 dV _ H 3 / 2 
~dt + 3Hd^ ~ 2tt ^ 



(4) 



where the noise field £ is defined in such a way that its 
correlation function simply reads 



(5) 



where 5{z) is the Dirac function. The normalization 
of the correlation function is chosen in order to re- 
produce, for a free field, the ordinary de Sitter result 
( lp 2 )=HH/(4n 2 ). 

At this point, two situations are possible, leading to 
very different technical problems. The first possibility 
corresponds to the case where the scalar field is a test 
field in a fixed background. This means that the factors 
H which appear into the Langevin equation (at the de- 
nominator in the second term and at the numerator in 
the third term) must be considered as functions of time 
but not as functions of the coarse-grained field. This is 
obviously an important simplification and, in this case, 
the noise is said to be non-multiplicative. In such a situ- 
ation, the derivation of the Langevin equation is unam- 
biguous and on a firm basis. In order to see how the 
formalism works, let us quickly consider the case where 
H is constant in time (i.e. de Sitter background) and 
V(<fi) = m 2 cj) 2 /2. Then, the solution of the Langevin 
equation can be found explicitly yielding 



cp(t) 



,-m 2 (t-i ln )/(3ff) 



#3/2 

2tt 



tin 



,(6) 



where (p- m is the initial value of the field. Then, we can 
easily deduce the two-point function and we obtain 



3H 4 
87r 2 m 2 



3H 4 
8n 2 m 2 



3 -2m 2 (t-t in )/(3/f) 



After a transitory regime, one sees that the two-point cor- 
relation function goes to (y> 2 (t)) — > 3iJ 4 /(87r 2 m 2 ). This 
well-known result has already been obtained, for instance 
in Ref. 0] by solving the Fokker-Planck equation. Let 
us notice that we could have also found the solution and 
computed the correlation function in the case of a colored 
noise. 

The second situation corresponds to the situation 
where one takes into account the back-reaction of the 
coarse-grained field on the geometry. Technically, this 
means that the Hubble parameter in Eq. (0J becomes 
a function of the coarse-grained field itself. In other 
words, the noise becomes multiplicative. In this case, 
it is necessary to have one more equation and one natu- 
rally assumes that the Friedman equation (in the slow- 
roll approximation) holds for the coarse-grained quanti- 
ties, namely 



3m' 



-V{ V ) = -V{v). 



(7) 



This case if of course the most interesting since it corre- 
sponds to the case of inflation. The coarse-grained field 
becomes the coarse-grained inflaton which drives the evo- 
lution of the background. 

Unfortunately, as is well-known in the case of a mul- 
tiplicative noise, the derivation of the Langevin equa- 
tion becomes also more problematic, see for instance 
Refs. [^j. Roughly speaking, this is due to the fol- 
lowing. When we promote the field <p to a stochastic 
quantity, there is some arbitrariness in defining the term 
i? 3 / 2 £ in the Langevin equation. Indeed, the quantity 
jj3/2 or jgi na tes from two contributions. The first one 
comes from the term H present in the damping term of 
the Klein-Gordon equation which, according to the rule 
outlined above, should be promoted to a stochastic quan- 
tity. The second one comes from the normalization of the 
noise correlation function which is an ordinary function. 
Therefore, the problem arises because one could choose 
to promote a different power of the Hubble constant to 
a stochastic quantity, say H 3 / 2 ~ x and put the remain- 
ing term, H x , into the normalization of the noise. In 
this case, the Langevin equation would lead to differ- 
ent results. In the present paper, we consider H 3 ^ 2 as a 
stochastic quantity. Finally, we notice that there is also 
the ambiguity in the choice of the calculus. Here, we 
work with the Stratonovitch calculus. 



III. INFLATION 
A. Classical Evolution 

Having specified what our basic set up is, we now turn 
to the application of this formalism to inflation. For sim- 
plicity, in the following, we restrict ourselves to single- 
field "chaotic" inflation models It turns out useful 
to parameterize the potential in term of the dimension- 
less scalar field tp/m pl . Explicitly, we take (with n > 2) 
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V[fp) = M — 



(9) 



The Cosmic Microwave Background Radiation (CMBR) 
anisotropy observations constrain the value of Vo- For 
small £, the multipole moments are given by 



C,= 



2H 2 



1 



25eml t £(£ ■ 



1) 



(10) 



and what has been actually measured by the COsmic 
Background Explorer (COBE) [22] and the Wilkinson 
Microwave Anisotropy Probe (WMAP) Q satellites is 
QLs-ps/T 2 = 5C 2 /(4tt) ~ 18 x 10- 6 /(2.7) .6x lO" 6 . 
The quantity e is the first slow-roll parameter [23] and for 
chaotic models, it reads e ~ n/{AN^, + n) where iV* ~ 60 
is the number of e-folds between Hubble radius exit and 
the end of inflation. Putting everything together, we find 
that Vq is given by 



m4 



90 



16 



(47V* + n)"/2+i ^ n 



n/2 



?rms-PS 



(11) 



From an observational point of view, all the models such 
that n > 5 are now excluded by the WMAP data, the 
quartic case being on the border line, see Ref. [24|. 

For the simple potentials considered here, the slow-roll 
equations can be integrated exactly. For this purpose, it 
is convenient to express everything in terms of the total 
number of e-folds defined by 



N = In | — 

. Clin 



(12) 



such that, initially, one has N — 0. Then, the classical 
field, i.e. the solution to the slow-roll equations of motion 
without the noise, reads 



<Pc\ 




4?r 



(13) 



where (p c \(N = 0) = Lp- m . The model remains under con- 
trol only if the energy density is below the Planck energy 
density. This amounts to the following constraint on the 
initial conditions ip in /m rl < (m't /Vo) 1 ' n . Inflation stops 
when the slow-roll parameter e is equal to unity corre- 
sponding to </2 cn d/m- p , = n/(4\/7r). As a consequence, 
one can easily check that the argument of the square 
root in Eq. (|13|l remains always positive. Finally, the to- 
tal number of e-folds during inflation is simply given by 
N T = 4ir(ipi n /m pl ) 2 /n — n/4. This number can be huge 
if the initial energy density of the inflaton field is close 
to the Planck energy density. 



Therefore, we use perturbative techniques. We consider 
the coarse-grained field tp as a perturbation of the solu- 
tion ip c i of the classical equation, i.e. we write 



<p(t) = <p c i(t) + 6<pi(t) + 5<p 2 (t) 



(14) 



where the term Stpi (t) depends on the noise at the power 
i. Clearly, this expansion is valid as long as Sif2 < Sipi < 
tp c \. Expanding up to second order in the equation of 
motion, we get two linear differential equations for Stpi 
and 8if2, namely 



and 



— — + -H (tpcijoipi = 

at k 2ir 



— — + -H {(pcl)0tp2 = 0<Pi 

at k k 



(15) 



4tt 



(16) 



where a prime denotes a derivative with respect to the 
field. These equations can be solved by varying the 
integration constant. Let us first consider the equa- 
tion for 8ip\. If the initial conditions are such that 
6ipi(t = tj n ) = 0, then the solution reads 



6<px(t)= /dr g3/2(yicl) exp[- [da 

J tin 2lT L Jt 



_2H"(<p cl ) 



This expression can be further simplified. If we use the 
classical equation of motion, then one can write the ex- 
ponential term as 



exp 



2H" 
da 



exp 



¥>cl(0 JJII 

dip lF 



¥>ci(t) 



H'Mt)} 



H> [<p( r )} ■ 
(18) 

Inserting the above expression into Eq. (|17fl . one finally 
arrives at 



2vr 



H' bci(r)] 



(19) 



We are now in a position where the various correlation 
functions can be calculated exactly. Since 8ipx is linear 
in the noise £, the mean value obviously vanishes 



(S(fi) =0. 



(20) 



Let us now evaluate the two-point correlation function 
calculated at the same time, i.e. the variance. Making 
use of Eq. JSJ, one obtains 



Perturbative Solutions 



-(-) 

2 \2n J 



fin 

dip 

Vol 



(21) 



In general, the Langevin equation cannot be solved an- 
alytically even for the simple potentials given by Eq. ^ . 



So far, the expressions presented above are general and 
do not rely on a particular shape of the inflaton potential. 
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FIG. 1: Evolution of (5ip\) /ip 2 ci versus the number of e-folds during inflation for three different chaotic potentials characterized 
by n = 2, 4, 6. The initial values of the inflaton field corresponds to an initial energy density of V ln = 0.5rn Fl and implies that 
the total number of e-folds during inflation is not the same for the three models under considerations. This is also the reason 
why the curves do not stop at the same place. For the n = 2 model, (5<pi) always increases and, at some point, the perturbative 
scheme used to calculate it breaks down. For n = 4, (8<p1) /ip^ tends to a constant that is difficult to visualize because of the 
e-folds logarithmic scale. For n > 6, {Sipf) possesses a maximum during inflation and then decreases. 



If we now specify the calculation to the chaotic potential 
given by Eq. then the variance takes the form 



4 V(ip cl ) 
3n m 4 . 



m , 



m , 



fir 
?7l 



(22) 

Since ip c \ is always smaller than ip\ u (because the field 
rolls down its potential) , the above quantity is increasing 
with time and always positive as required. 

We now turn to the equation of motion for the second 
order perturbation 8(f2- It can be solved by following 
exactly the steps that were described before. Then, the 
solution can be written as 



TTl ft TTlll 



3H' 

47T 



(23) 



As expected the second order perturbation is quadratic in 
the noise. One can now easily evaluate the mean value of 
Sif2{t), taking into account a factor 1/2 which originates 
from the fact that the Dirac ^-function appearing in the 
noise correlation function is centered on an integration 
limit. One obtains 



(Sip 2 ) n + 2 V(tfcO ( <A;i \ ™ /2 1 



3n m 4 

p; 

n-2 V(<Pc\) f <£h 
3n to 4 V 



m pi 

-n/2-l 



i/2+2 



gel 



m F 
n/2-2 



t/2+2 



n/2-2' 



(24) 
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The quantities (<Vi) /Pel anc ^ (^2) /Pel are represented 
in Fig. and |2 for different potentials, i.e. for different 
values of the power index n. We find that for n > 4 both 
quantities reach a maximum after an increasing phase, 
and then decrease to zero at late times, always remaining 
smaller than unity for any choice of the initial conditions 
(thus fully validating the perturbative treatment). Con- 
versely, we see that their late time behavior is constant 
for n = 4 or rapidly increasing for n < 4, and appears 
in both cases to violate the validity of the perturbative 
approach for the extreme choice of an initial condition 
corresponding to an energy density near the Planck scale. 



C. Classicalization 

Let us study the late time behavior of (61P2), as given 
by Eq. (|24p. in more details. In this regime, one has 
ifd <C Pin- Then, it follows that 



(Sip 2 ) ( Pel 

oc 

m n , \ m„, 



i/2-l 



(25) 



From the above equation, one deduces that we have ex- 
act classicalization for the cases n = 2 or n = 4. Indeed, 
in this case, we have (p) ~ p c i + (S1P2) oc p c i which is, by 
definition, what we mean by classicalization. For n > 4, 
the situation is slightly different. In this case, (p) is not 
exactly proportional to ip c \ but the extra contribution 
(S(f2) decays faster than ip c \ and, hence, becomes negligi- 
ble as can also be checked in Fig. [3 In this way, we also 
recover a classical behavior. 

Let us now see what happens to a trajectory which is 
far from the mean value. From Eq. I|22[l. it is easy to 
check that, at late times, we also have \J (<5p 2 ) /m pl oc 
(<y9 c i/m pl )"/ 2_1 (a quantity which gives an estimate of 
the amplitude of the probability distribution). It follows 
than not only the mean value of the distribution but also 
its width evolves according to the classical solution. 

Finally, the following remark is in order. The fact that 
the quantum effects are negligible does not mean that, 
at late times, the value of the quintessence field (ip) is 
the same as the corresponding classical field value calcu- 
lated with the same initial conditions. Since the quantum 
could have accumulated in the early phase of evolution 
causing a deviation from the classical trajectory, this ac- 
tually means that the late time value of (ip) can only be 
viewed as p c i but calculated with different initial condi- 
tions. 



D. Comparison with Other Methods 

In this section, we compare the results of the previous 
section with some known results of the literature. 

In Refs. [2M [2(j, it has been shown that there is a 
case where the Langevin equation can be solved exactly. 
In order to permit a more direct comparison with those 



studies, we write the constant Vb which appears in the 
expression of the potential © as Vo/m* — 3A„/(87r), 
i.e. in terms of the coupling constant A„. Then, the rea- 
soning goes as follows. The Langevin equation can be 
written as 



A (JL 

dt \m Pl 
,3/4 

An 



2tt 



,1/2 



3n/4 



1/2-1 



_£_ 



(26) 



There is an exact solution if n = 4 because, in this case, 
the Langevin equation takes the form of a Bernoulli equa- 
tion. The simplification comes from the fact that the 
equation can be put under the form of a linear equation 
for some power of the field (in the present context, this 
power is —2 hence the form of the next equation). Let 
us also notice that the choice n = 4/3 would also lead 
to a Bernoulli equation but this case looks less interest- 
ing since the power index of the potential is no longer 
an integer. The Bernoulli equation can be solved exactly 
because it may be brought into a linear form by a change 
of variable. The solution can be written as 



-2 

V \ _ e \AIro p! (t-t to )/7r 
m, 



pi 



771„ 



-2 A 3/4 



4 m 1 / 2 



< I dre 



-\A4TOp, (T-tuO/TT 



(27) 



Taking the inverse square root of the general solution one 
then obtains 



ip(t)=ip d (t) [1-*(<)P 1/2 , 



(28) 



where ip c i(t) is the classical solution that can be expressed 
explicitly in terms of cosmic time 



<Pci(t) = Pin exp 



2n 



(t - tin) 



(29) 



while the quantity *5f(t) is defined by 



*(*) 



,3/4 

A 4 1/2 

— — raj 





>cl(r)" 







£(r) . (30) 



This function can be treated as a new dimensionless 
Gaussian noise with vanishing mean value and whose 
variance reads 



A4 

2tt 



gcl 



(31) 



We are now in a position where the various correlation 
functions of the scalar field can be computed. Here, 
we compute the two-point correlation function. Using 
Eq. QBJ , one gets 




Pel 



2 +00 



(32) 



fe=0 
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FIG. 2: Evolution of (8ip2) /(fi c \ with the number of e-folds during inflation for three different potentials. The initial values of 
the scalar field always correspond to an initial energy density of Vi n = O^m^ . As a consequence, the total number of e-folds 
during inflation is not the same in these three models and this explains why the three curves stop at different N. A value of 
(6<p2) /Vci greater than one signals a breakdown of the perturbative approach. For instance, this is the case for the n = 2 model 
at the end of inflation. 



Let us notice that we could have also computed any 
(i.e. n-point) correlation function of the field (including of 
course the mean value) using the same technique. Since 
we have to deal with a Gaussian process, only the even 
correlation function are non vanishing. In addition, for 
a Gaussian process we have ^ 2k ^ = [2k — 1)!! (\I/ 2 ) fe . 
Therefore, we obtain a series the general term of which 
can be written as 2 k T (k + 1/2) (^ 2 ) k A/7T, where we 
have used Eq. (8.339.2) of Ref. [23. At this point two re- 
marks are in order. Firstly, it is easy to convince oneself 
that the series (|32|) is in fact divergent. The interpre- 
tation of this fact is of course a little bit problematic 
but, on the other hand, it is well-known that perfectly 
well-defined distribution functions can have no moments. 
Below, we also give another interpretation of this fact. 
Secondly, looking at the expression of the new Gaussian 
noise *]/, we see that the series is in fact an expansion in 
the coupling constant A4. Since this one is in fact tiny, 
one can try to work with the first term of the series only. 
Then, one obtains 



<gcl 



(33) 



This expression is in fact exactly similar to the one ob- 
tained previously for (p c \(t)+S(pi(t)+S(f2{t) and therefore 
leads to the same correlation functions as before. It is 
easy to show that the expansion in the coupling constant 
is nothing but the expansion in the noise used before. 



Our method is more general because it is not restricted 
to the case n = 4. This is because the expansion is di- 
rectly performed in the Langevin equation rather than in 
its solution. The drawback of this last method is clearly 
that it is first necessary to find a solution of the Langevin 
equation, which is not an easy task, before the expansion 
can be taken. Let us emphasize again that, even if an 
exact solution is known, an expansion is still necessary 
because only some power of the field is generally obtained 
(in the present case <p~ 2 ) and, in order to find the expres- 
sion of the stochastic field itself, one should then compute 
the root of the solution. 

The same conclusions can be obtained by means of a 
model where the random walk of the field is modified 
by the presence of a reflecting barrier. This model is 
explored in the Appendix lAl 

Another possibility studied in the literature is the use 
of the so-called scaling solutions, see Ref. . Let 
us quickly review the method. The first step con- 
sists in rendering the stochastic process non multiplica- 
tive. This can be done by means of the transforma- 
tion x = 16 x 3 3/4 tt 4 V- 3/4 (6)d8 which reduces the 

Langevin equation to dx/dt = — dV/dx + £(t) where 
V = -3x(4TT 2 ) 5 /V[tp(x)]. Explicitly, one has -dV/dx = 
g x (ie/to p ,)(" +4 )/( 3 "~ 4 ) where g is a constant which de- 
pends on Vo and n. The next step is to consider the 
time-dependent nonlinear transformation 



t]( x ,t) = F^ie-^-^Fix)} , 



(34) 
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where the function F is defined by, see Ref. [28 

v 2(n-4)/(3n-4) 

L_ I 'IV ' 1 ~ ^) I X ' 

m„, I 2o(n - 4) V 



2g(n - 4) 

2(n-4)/(3n— 4)-i 



(35) 



with 7 = .g(x in /m pi ) 2 ( 4 -™)/( 3Tl - 4 ). Then, it is straight- 
forward to show that the new stochastic process r\ obeys 
the following equation 



dij 
dt 



n — 4 
l + 2ff- At -tin) 



XT) 



3n - 4 



2(4-n)/(3ra-4) 



-(n+4)/[2(n-4)] 



£(t). (36) 



So far, everything is exact. However, the above Langevin 
equation cannot be solved exactly. The so-called scaling 
solutions correspond to replacing r) by r) c \ in the right 
hand side of Eq. (|36|l . Then, clearly, the Langevin equa- 
tion can be integrated in this approximation. Therefore, 
the scaling solutions are nothing but the result of an ex- 
pansion in the noise and, again, are similar to the solu- 
tions found at the beginning of this section. We notice, 
as shown in Ref. |2S| . that the quantity rj c \ is in fact a 
constant. This is consistent with the fact that this is 
the zeroth order solution (i.e. for which the noise term 
is simply neglected) of the above-mentioned expansion 
for which the Langevin equation simply reads dr)/dt = 0. 
Finally, in Ref. 28] , a saddle point approximation is used 
in order to calculate the effective dispersion. The result 
found is, see Eq. (28) of Ref. 

(ip 2 ) cS cx ff 2 (^ cl )sinh[lntf- 4 /"(^ cl )] (37) 
= p« (<p 2 + ^ 2 ) , (38) 

which is exactly the result obtained in Eq. I|22l) . This 
reinforces the result that the scaling solutions are very 
similar or even identical to the solutions exhibited here 
by means of the expansion in the noise term. 

In conclusion, the solutions obtained previously for the 
stochastic innaton are explicit and consistent with those 
already found in the literature. In the sequel, we use 
them as a description of the background in which the 
quintessence field lives. 



IV. QUINTESSENCE 

A. Classical Evolution 

In order to explain the accelerated expansion of the 
universe, one postulates the presence of the quintessence 
field Q. This field is a test field during almost all the cos- 
mic evolution and becomes dominant only recently when 
it drives the accelerated expansion. As is well-known 



for a scalar field, the equation of state ujq — Pq/pq is 
time-dependent and, crucially, can be negative. The de- 
tailed evolution of luq clearly depends on the shape of 
the quintessence potential W(Q). An interesting choice 
is the inverse power law potential (with a > 0) which 
was first studied by Ratra and Peebles in Ref. 



W(Q) = Wo [ — 



(39) 



During the radiation dominated era, it is possible to 
find an exact solution of the corresponding Klein Gordon 
equation for which Q oc a 4 /( Q+2 ) or pq oc a^ 4a ^ a+2 \ 
This is also possible for the matter dominated era for 



which one has Q oc a 3 /( Q+2 ) or 



Pq oc a 



-3a/(a+2)^ Thg 



two solutions we have just mentioned can also be ex- 
pressed by means of the following equation p| 



d 2 W(Q) 
dQ 2 



9a 
2 _ 



1 



(1 



(40) 



this expression being valid both during the radiation and 
matter dominated epochs. We can re-write the param- 
eter ujq as ujq = (auj B — 2)/(a + 2) where lo b is the 
equation of state of the background, i.e. either 1/3 or 
0. Since pq redshifts slower than the background energy 
density, the scalar field contribution will eventually be- 
come dominant. From the above equation, it is easy to 
see that, when the quintessence field is about to domi- 
nate, its value is in fact of the order of the Planck mass. 
Then, the value of Wq is constrained by the fact that the 
quintessence energy density is almost the critical energy 
density today. This gives 



Wp 



8ir to 2 



(41) 



where Hq is the Hubble parameter today, i.e. Hq 



-61, 



Equipped with this relation, one can also es- 
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timate the ratio Wq/Vq and one gets 



Wo 
V 



(4A„ 



W2+1 



240tt 



(—) 

view 



«/2 H, 



Qr 



ms-PS 

T 



This number is obviously extremely small, ~ 0(10 no ). 

The main property of the solutions described before is 
that they are attractors @- This means that there is no 
need to fine-tune the initial conditions and that the so- 
lution will be on tracks today for a large range of initial 
conditions. Let us be more precise about this particular 
point. Usually, one fixes the initial conditions at the end 
of inflation, at a rcdshift of z = 10 28 . Then, the allowed 
initial values for the energy density are approximatively 
such that H El EH: 10" 37 GeV 4 < p Q < 10 61 GeV 4 
where 10 _37 GeV 4 is the background energy density at 
equality whereas 10 61 GeV 4 represents the background 
energy density (i.e. the radiation energy density) at the 
initial redshift. If, for instance, a — 6 and if the scalar 
field starts at rest, this means that the initial values of 
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the field are such that 10 _18 m pi < Q m < 10~ 2 m Pi just 
after inflation. If Q ln is large initially, then the attractor 
is joined quite recently. Of course, the range of allowed 
initial conditions, when expressed in terms of the initial 
field, depends on the value of the parameter a. 

The values of a are constrained by the measurement 
of the equation of state today. The larger a is, the larger 
the equation of state parameter uq is today. Since it 
is known that ljq cannot be too different from —1, this 
means that a cannot be too large. In fact, this conclu- 
sion rests on the particular shape of the Ratra-Peebles 
potential. From a model building point of view, it is 
more natural to consider the SUGRA potential given 

by OHHE0 



W(Q) = W e 



Q 



(43) 



Then, the equation of state parameter today is modified. 
Since, today, the value of the field is of the order of the 
Planck mass, the supergravity exponential factor plays 
an important role in this regime. This has two conse- 
quences. Firstly, the equation of state is pushed toward 
the value — 1 because the exponential factor increases the 
importance of the potential energy with respect to the 
kinetic energy. Secondly, the value of loq becomes al- 
most independent of a because, again, the exponential 
factor dominates. As a consequence the constraint on a 
can be relaxed and any value is in fact a priori possible. 
Moreover, in the very early universe, one has Q <C m pi 
and, this time, the exponential factor becomes one. In 
this case, the SUGRA potential has the same shape as 
the Ratra-Peebles potential. In summary, if we have the 
SUGRA model in mind (which is, from a high-energy 
point of view, well- motivated), then, in the early uni- 
verse, we can safely work with the Ratra-Peebles as an 
effective model (which is simpler, technically speaking) 
but without the usual restrictions on the parameter a. 

We have seen that the initial conditions are usually 
fixed just after inflation, at the beginning of the radiation 
dominated era. In this paper, we study the behavior of 
the quintessence field during the phase of inflation itself. 
One of our main goal is to check whether the behavior 
of Q during inflation is compatible with the allowed ini- 
tial conditions (in fact, "final" conditions from the point 
of view of the present study) described before. Clearly, 
if the final value of Q (at the end of inflation) is not in 
the allowed range, then the quintessential scenario is in 
trouble. In fact, in order to avoid the above-mentioned 
situation, our study rather helps us to put constraints on 
the initial conditions, not at the beginning of the radia- 
tion dominated era as before, but at the beginning of the 
inflationary phase. 

We start with a study of the classical evolution of the 
quintessence field (i.e. without the quantum effects). The 
classical quintessence field obeys the usual equation of 
motion for a scalar field in a FLRW background with 
the Hubble parameter depending on the classical inflaton 



field </? c i, namely 



Q cl + 3H(^ cl )Q cl + W'{Q cl ) = 0. 



(44) 



Whenever the friction term is large, one can neglect the 
double derivative term, obtaining the standard slow-roll 
equation 



dQci W'jQci) 
dt 3H(<p cl ) 



= 0. 



(45) 



The consistency of this assumption can be directly 
checked. Taking the time derivative of the last equation 
one can actually show that 



Qc 



H((f c \)Q c 



W"(Q d ) 



(46) 



where e = —H/H 2 is the usual inflaton slow-roll param- 
eter which is by assumption much smaller than 1 . There- 
fore, the slow roll approximation can be applied to the 
equations describing the motion of the quintessence field 
whenever the following condition is satisfied 



W"(Q cl ) 
3# 2 (¥>ci) 



«1, 



(47) 



If one applies this condition to the Ratra-Peebles poten- 
tial, one arrives at 



Qd 



a+2 



» 



a{a + 1)W Q 
SttV( Vc1 ) 



(48) 



Let us now turn to the solution of Eq. I|45|l . For the 
Ratra-Peebles potential given by Eq. i|39|) . the solution 
to this slow-roll equation can be expressed as 



Oh 



a+2 



ot(a + 2)W( 



l/{a+2) 



(49) 



Using the chaotic inflationary potential © we get for any 
value of the power index such that n > 2 




a+2 



X 



a(a + 2) W 
n(2-n) ~Vo 

2-n / \ 2-n' 



l/(a+2) 



(50) 



while, for n — 2, the result reads 

a+2 



Qci 



Oh 
rn r 



a(a + 2) W , if 



In 



Vo ificl(t) 



l/(a+2) 



(51) 

These results are consistent with those obtained in 
Ref. The evolution of the quintessence field in the 

slow-roll regime for different initial values is plotted in 
Fig. |3 Let us quickly discuss these solutions. One can 
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FIG. 3: Classical evolution of the quintessence field for different initial conditions. A massive chaotic inflaton potential 
(i.e. n — 2) has been assumed and the parameter a in the Ratra-Peebles potential has been taken to be equal to 6. The initial 
condition for the inflaton field is such that ip- ln = 2.5 x 10 6 m pi corresponding to an initial energy density of Vi n = O.lm^ . The 
green dotted-dashed line signals the limit of the slow-roll approximation. Below this line, the slow-roll approximation is no 
longer valid. 



obtain a more compact expression if one notices that 
(fd ~ ^m[l — N/(2N T )] where N T is the total number 
of e-folds during inflation. This expression is valid pro- 
vided that N < An(ip iri /m pl ) 2 /n — N T + n/4 ~ N T . 
Then, one has 



a(a + 2) W ( <p 
2nV 



/ Qin} 


-a-2 N - 


l/(a+2 









(52) 



If the initial value is large enough, then the dynamical 
term will remain negligible and the field is frozen until 
inflation ends. As a matter of fact, for any value of n we 
obtain that Q c \{t) — Qin at all times if the initial value 
satisfies the constraint 

a+2 



III, 



> 



Wo 
V 



(53) 



If the above condition is satisfied at initial time, then, ob- 
viously, the condition given by Eq. (|48|l is also satisfied. 
There is also an intermediate regime for which the condi- 
tion (|53|l is not satisfied but (|48|l is. In this case, the field 
is not frozen, even at initial times. Finally, if Q ln is suffi- 
ciently small so that the condition (|48|l is violated, then 
we expect the quintessence dynamics to rapidly bring 
back the field to the slow-roll regime, where the previous 
considerations apply. 



B. Perturbative Solutions 

As mentioned at the beginning of this article, in or- 
der to study the evolution of the quintessence field, it is 
not sufficient to integrate the classical equation of mo- 
tion since the quantum effects can play an important 
role and modify the classical evolution. We now ana- 
lyze the stochastic behavior of the quintessence field, in 
the case where the total energy density is still dominated 
by the vacuum energy of the inflaton field. The stochas- 
tic evolution of the quintessence field Q is controlled by a 
Langevin equation which, in the slow-roll approximation, 
reads 



dQ 
dt 



W'{Q) H 3 / 2 (v) 



2tt 



(54) 



3H(<p) 

where £q is another white-noise field such that 

mm Q (t')) = s(t-t'), <£ Q (W)>=o. (55) 

The solution of the Langevin equation l|54|) depends ex- 
plicitly on £q but also on the inflaton noise £ through the 
coarse-grained field (p. 

In order to find an approximate solution to Eq. 154fl . 
one may try to use the same perturbative technique as 
the one used before for the inflaton field. Therefore, we 
expand the quintessence field about the classical slow-roll 
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solution (|49|l and write for the perturbed quantities 8Q\ and SQ 2 are given by 

the following expressions 

Q(t) = Qd(*) + 6Q1 + 8Q2 + ■ ■ ■ • (56) 
Then, it is easy to establish that the equations of motion 
I 

dSQ 1 W^q^) W'{Q ci )H'{ Vci ) x ff 3/2 (y> cl ) , (vr . 

— + im^) 5Ql = m^) Vl + ' (57) 

dSQ 2 , W"(Q cl ), n W'(Q d )H'(<p cl ) s W'(Q cl )H"(^ cl ) 2 W'(Q cl )H> 2 (<p cl ) x 2 



dt 3H{tp d ) " z 3H 2 (<p cl ) rz 6H 2 (<p cl ) rL 3ff 3 (<A 



+ 3ff%a) ^ lgQl " Ql + 4^ 6 ^ Q ■ (58) 



Although these equations look quite complicated, they 
can be solved easily because (by definition) they are lin- 
ear. The solution for 8Q\ reads 



SQi(t) = W'(Q C 



tin 



3ff 2 (^ cl ) 

. g 3/2 (^d) 
2ttW'(Q c1 ) 



8<Pi(t) 



(59) 



and, as required, is linear both in the quintessence noise 
£q and (through Stpi) in the inflaton noise £. As a con- 
sequence, 5Qi has a vanishing mean value 



(SQ 1 )=0, 



(60) 



but a non- vanishing variance given by the sum of two con- 
tributions originating from the inflaton and quintessence 
noise variances, namely 



(SQi) = 



W' 2 (Q cl ) 
9 

, W' 2 (Q cl ) 
in 2 



tin Jt 



H'{t) H'( V ) 
, b H 2 {t) H 2 (r]) 
1 '* H 3 / 2 (r) H 3 / 2 ( V ) 



W'{t) W'{rj) 



{5tp\{r)8(pi{ri)) drAr] 

^Q(r)^ Q (v))dTd V 



= (SQi) b + (SQi) k Q 



(61) 
(62) 



Let us notice that there is no mixed contribution since the cross-correlation (££q) = 0. The detailed calculation of 
(SQ 2 ), in particular its explicit expression in terms of the inflaton field and/or the number of e-folds N, is rather 
lengthy and is carried out in the Appendix FBI 

Let us now turn to the second order correction. The solution for 8Q 2 can be written as 



SQ 2 (t) 



W'(Q c i) 

W"'(Q ci ) 



/ ff'teci) 



W(Q A )H(tp 



SQi(r 



5ip 2 {r 



H"{<p cl ) H' 2 (<p cl ) 



6H 2 (ip cl ) 3H 3 (<p cl ) 
3 H l / 2 {y cl )H'{y c {) 



s 2r ^ , W"{Q cX )H'{<p cl ) x 
SViiV + ^ TT ■ ? , — sSvi{t)8Qx(t) 



3W'(Q ci )H 2 ^ cl ) 



■in 



W'(Q C {) 



S<Pi(t)£q >dr. 



(63) 



As expected, one sees that SQ 2 is quadratic in the noises. From the above expression, one deduces that the mean 
value of SQ 2 is non-vanishing and is the sum of various terms 



(6Q 2 ) = (8Q 2 ) \s V2 + (SQ 2 ) \srf + (SQ 2 ) \s VlSQl + (8Q 2 ) \sqk^ 

I 



{8Q2 



(64) 



where the last term in Eq. I|63|) does not contribute be- cause 



(65) 
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If we had not taken into account the stochastic behavior 
of the inflaton, only the term (8Q2) l<5Q 2 (£ Q ) would have 
contributed. Again, the explicit expressions of each term 
are given in the Appendix iBl 

Let us now quickly present what the outcome of the 
perturbative approach applied to the evolution of the 
quintessence field is. The main result is that when the 
classical evolution of the quintessence field is negligible, 
i.e. when the condition (|53|l holds, the variance reads 



16m 



Vn 



3n(n + 2) mi 



/ \ n+2 




( ) 









n+2 



1-1 (66) 

which is the same result already obtained in Ref. [l5j that 
we recover here with a different method. 

Another general result that is established in details in 
the Appendix IBl is that, when the perturbative approach 
is valid, the contribution coming from the inflaton noise 
is completely negligible, see for instance Figs. |S1 

Unfortunately, this is not a very strong argument be- 
cause it is possible to see that the perturbative treatment 
is well under control only in a very small region of the 
parameter space. Indeed, one must check that the con- 
ditions 



(SQf) 



< 1. 



(SQ 2 ) 
Qci 



< l. 



(67) 



are fulfilled if one wants the perturbative approach to be 
valid. This leads to several constraints that are summa- 
rized in Fig. More precisely, from Eq. I|62l) , one sees 
that the first condition leads to two constraints while, 
looking at Eq. Q64p. one notices that the second one leads 
to fives constraints. As is apparent from Fig.0J the most 
stringent constraint comes from the fact that the variance 
originating from the quintessence noise must be small in 
comparison with Q 2 V Explicitly, working out this condi- 
tion, it boils down to 



Or 



> 



(fir. 



n/2+1 



(68) 



Let us also notice that this condition is a much more 
stringent condition than (|53J) . see Fig. 01 For most val- 
ues of Q[ n and ip- m , we rapidly get that (8Q\) jQ\ x 1, 
meaning that a large part of the initial condition space 
cannot be described by means of the perturbative ap- 
proach. Therefore, there is the need for a different ap- 
proach. 



C. The Reflecting Wall 

The failure of the perturbative treatment is a signal of 
the fact that the classical evolution is not a good zeroth 
order solution for the perturbative expansion. The rea- 
son is that the equation of motion is dominated by the 
diffusive term due to the noise while the classical drift is 



in fact sub-dominant. It is therefore natural to neglect 
the classical term and to solve the approximate equation 



dQ 
dt 



2tt 



(69) 



The solution to this equation can be written as 

Q = Qin + I dr n _ W £(r) • (70) 



2tt 



As a first step, we also neglect the inflaton fluctuations 
and take ip to be the classical field tp c \. As a consequence, 
Q has a Gaussian probability distribution the mean of 
which is given by (Q) = Q m with a variance <Jq = (Q 2 ) — 

(Q) 2 which can be expressed as 



2/+\ — 



dr^ 



4tt 2 



16m 2 



Vo 



3n(n + 2) m 4 



/ \ n+2 




( ) 


{ \ 







n+2 



.(71) 



One recognizes Eq. (|66|l . first obtained in Ref. and 
that we have re-derived in the preceding subsection by 
means of the perturbative approach. Here, we have 
shown that Eq. (|66|l can be valid even when the per- 
turbative approach breaks down. 

At this point a remark on the notation is in order. 
In the following Oq always denotes the variance of the 
quintessence field while a 1 is just a function defined by 
the above expression which turns out to be equal to <Tq 
in the situation where the classical drift and the inflaton 
noise are neglected. 

The above model is in fact too simple for the follow- 
ing reason. The function <7o(£) increases with time and, 
at some point, is larger than Q ln . In this case, there is 
a finite probability (tending to 50% at late times) that 
the quintessence field becomes negative. Clearly, this is 
not possible because, in this case, the classical term in 
the Langevin equation becomes dominant and prevents 
Q to become negative. In other words, in this regime, 
the classical drift cannot be neglected. So, it seems that 
we are in fact back to the original problem which consists 
in solving exactly the full Langevin equation. However, 
there is a simple way out. Indeed, we can model the ef- 
fect of the classical term by considering that there is a 
perfectly reflecting wall at Q = the purpose of which 
is of course to prevent the quintessence field to become 
negative. As discussed in Ref. [2j|, the probability dis- 
tribution of a random walk is modified by a reflecting 
barrier in a way that is easy to estimate. Let us assume 
that we start with a normalized probability distribution 
for the variable x, P(x)dx — 1. Then, let us put a 
reflecting wall at x — a such that only the values x > a 
are allowed. Then, Ref. 29] tells us that the new proba- 
bility distribution is P(x) + P(2a — x) and it is easy to 
check by mean of a simple change of variable that it is 
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FIG. 4: Constraints coming from the requirement that the perturbative approach is valid for the case n = 2 and a = 6. 
The values of the initial conditions y>i n and Qi n are chosen such that the corresponding energy density is never Planckian, 
i.e. V(<fii n ) < rn Pl and V(Qin) < m p ( ■ Each line corresponds to a condition specified in the figure and that we describe 
below. By convention, the allowed region is always the region above the corresponding line. The solid black line indicates 
the regime where Q c \ can be considered as constant, see Eq. |S3I- The dot-dashed green line corresponds to the region where 
(SQi) /Qci < 1) originating from the quintessence noise, is small, see Eq. 16811 . The dotted blue line corresponds to the same 
condition except that the variance now comes from the inflaton noise. Finally, the dashed red line corresponds to the condition 
{SQ2) /Qcl < 1 where the mean value originates from the quintessence noise. As we have seen before, see Eq. (1641 . there are four 
other contributions but one can easily show that they do not lead to new constraints and this is the reason why they are not 
represented in this figure. As discussed in the text, the most stringent constraint comes from the requirement (SQf) /Q% < 1. 



indeed normalized, i.e. j^°°[P(x) + P(2a — x)]Ax = 1. 
In the present context, we have a — and, hence, the 
resulting probability distribution becomes, for Q > 0, 



P(Q,t) 



-(Q- 



2 /(2-o 2 



-(Q + Q,n) 2 /(2<T ) 



27T<7 (t) 



(72) 



where ag is given by Eq. I|71|l which means that, for the 
moment, the contribution coming from the inflaton noise 
is still neglected. 

As mentioned before, the advantage of the reflecting 
wall model is that it permits simple analytical estimates 
of the relevant physical quantities. However, there is one 
feature of the model that is worth stressing here. The 
classical drift term which prevents the field to become 
negative depends on the parameter a but the wall, which 
is supposed to model this term, don't. Therefore, one 
limitation of the reflecting wall model is that we have 
lost the a-dependence of the result. Concretely, in the 
following, we will see that the mean value and/or the 
variance of Q are a independent. Only an exact solution 
(or a numerical calculation) could allow us to test the 
accuracy of this assumption. 

With the help of this probability distribution we 



can now calculate the mean and the variance of the 
quintessence field. For the mean, one obtains the fol- 
lowing analytical expression 



(Q) = to 



Qm Erf ( Q 
CO 



\%/2ct 



(73) 



where Erf (2) is the error function defined by Erf(z) = 
(2/a/7t) J* d£e~* . To our knowledge, this explicit for- 
mula is new and has not been given elsewhere. In the 
same manner, one has (Q 2 ) = + Qf n and, therefore, 
one has 



Q 



(QY 



(74) 



where (Q) is given by Eq. (|73J) . One notices the vari- 
ance of the quintessence field Oq is now different from 
the function CTg, Oq ^ (Tq. Let us emphasize again that 
the above results are not subject to the limitations of the 
perturbative approach. In particular, the difference be- 
tween the classical value and the quantum (or stochastic) 
average needs not to be small. 

If we insert Eq. JTlJ into Eqs. (JJSJ and CU, one ob- 
tains the mean value (Q) and the variance Oq as functions 
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FIG. 5: Evolution of (Q) (solid black line) calculated from Eq. 1731 (but with the inflaton noise taken into account) and 
compared with the corresponding classical evolution (dotted-dashed green line) Q c \. The dashed area represents the zone ±cr 
around the mean value where a has been obtained with the help of Eq. I|8ip. The difference between (Q) and Q c \ at the end of 
inflation can be many orders of magnitude (and is, in the present case, of about 8 orders of magnitude) demonstrating that it is 
crucial to take into account the quantum effects during inflation. The curves have been obtained for a massive chaotic model, 
i.e. n = 2 and for the Ratra-Peebles potential with a = 6. The initial value of the inflaton is ip- ln = 5.4 x 10 6 m pi corresponding 



to an initial energy density of VS n = 0.5 x m.p. 
Win = 1.25 x 10 _89 mi. 



The initial value of the quintessence field is Qin = 10 



corresponding to 



of the inflaton and/or the number of e- folds. At initial 
time, one has cto — » and therefore (Q) — > Q ln , where we 
have used the fact that Erf(z) — > 1 when z — > +oo. For 
the variance, one has in the same regime (i.e. initially for 

ifi -> (f m and (7 (*) < Qin) 

o$->og(t), (75) 

where we have used Eq. (I74|l and the behavior of the 
mean at initial times. This means that, in this regime, 
we recover the previous results obtained by means of the 
perturbative approach. 

One can also study what happens at late times when 
oo(t) S> Qin- In this case, one has 




The situation is reminiscent to the behavior of 
quintessence in the more traditional situation where the 
background evolution is dominated either by radiation 
or matter. Indeed, as it is the case in this well-studied 
context and as explained before, the late time evolution 
of (Q) is independent on the initial conditions, i.e. on 
Qin- In other words, there is an attractor for (Q) given 
by Eq. 1)76(1. In this situation, the final value of the 
quintessence field only depends on the initial value of 



the inflaton field (and, of course, on which kind of poten- 
tial is responsible for inflation: in the present context, 
it only depends on n). However, this conclusion should 
be toned down because we will see in the following that 
another quantity of interest, namely the probability that 
the quintessence field be on tracks today (or, equivalently, 
that its value at the end of inflation be in a given range) 
does depend on Q; n (and, in fact, also on <^i n ). 

Using the above expression, the variance of the 
quintessence field at late times can also be estimated. 
It is given by 

<4->( 1 -£)°3(*)- (77) 

We see that, as ao(t) increases and as the probability 
that the random-walking field is reflected by the wall 
becomes non negligible, the mean value (Q) is shifted 
toward higher values while its variance slightly shrinks. 

We are now in a position where we can come back to 
one of our starting questions, namely the influence of the 
inflaton noise. In particular, we study whether including 
the inflaton fluctuations can cause significant deviations 
from the above results and, if so, under which physical 
conditions this is the case. In order to take into account 
this effect, we re-start from Eq. {7UJ). In this formula, 
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FIG. 6: The quantity (Q) at the end of inflation versus the initial value of the inflaton field calculated according to Eq. 18311 
for two chaotic models, n — 2 and n = 4. As noticed in the text, (Q) at the end of inflation does not depend on the initial 
value of the quintessence field because there is an attractor solution given by Eq. I|76fl . The fact that the final value does not 
depend on the parameter a is due to the description of the classical drift term by means of the reflecting wall. The two curves 
start and stop at different values of the initial inflaton field because the minimal value (defined to be the one which leads to 
at least 60 e-folds) and the maximal value (defined to be the one which corresponds to V\ n = mi, ) are not the same according 
to the model (i.e. n) under considerations (for instance, </?i n = 3.1m p , for n — 2 while <p- m = 4.4m pi for n = 4). The horizontal 
lines represent the maximal values of (Q) (for different values of the index a) at the end of inflation such that the attractor is 
joined today. These values are computed by means of Eq. (I85H . 



the argument of the Hubble parameter is no longer the 
classical inflaton but the stochastic process studied in 
the previous section. Then, one can use the perturbative 
treatment studied before but, and this is of course the 
crucial point, only for the inflaton field. Indeed, we have 



seen before that this perturbative treatment (contrary to 
the perturbative treatment for the quintessence field) is 
almost always reliable. Performing a Taylor expansion of 
Eq. I|7L)|). one obtains 



(78) 



Equipped with this solution, one can now reintroduce 
the reflecting wall and recalculate the various mo- 
ments. Clearly, since the solution is still linear in 
the quintessence noise, the probability function of the 
quintessence field, taking into account the wall, is still 
given by Eq. I|72l) but, in the expression of P(Q,t), a 
should now be replaced by another function (of the clas- 
sical inflaton noise or of the number of e-folds) that, in 



the following, we simply denote by a (not to be confused 
with <jq). The replacement of <jq by a is the only change 
needed. Otherwise the expression of the new P(Q,t) is 
similar to the one given by Eq. I|72|) . Let us now deter- 
mine a explicitly. As before, a is simply given by the 
variance deduced from Eq. (|78)l using the fact that we 
have white noises. Eq. I|78|l implies that a 2 is given by 
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a 2 it) 



tin 



47T 2 



(79) 



that can now be easily integrated. The final result is 



where the mean value in the integral can be expressed 
through a power expansion of H 3 (ip) up to second order 
about its classical value yielding 



(H 3 { V )) = H 3 ( Vcl ) + 3H 2 H' {6<p 2 ) 
H{H'f + 



(8rf) , (80) 
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Let us now compare the function cr(N) with the func- 
tion 00(A). At early times, i.e. when ip c i(t) ~ ip- mi a 
linearization is sufficient in order to obtain a good ap- 
proximation (and this is in fact the case for a large part 
of inflation). A striking feature of the above result is that 
the extra term coming from the inflaton quantum fluc- 
tuations cancel out, the first non-vanishing contribution 
being of order 0[(N/N T ) 2 ], where we remind that iV T 
is the total number of e- folds. Therefore, for the main 
part of the inflationary era we can safely consider that 
a(t) ~ 0o (t). Again, we reach the conclusion that the 
effects originating from the inflaton noise do not play a 
crucial role. Conversely, at late times, when ip c \(t) <C (fi n , 
the above expression yields 



* 2 (t)=a 2 (t) 



1 



(n + 2)(7n + 4) V ir , 



3n(n + 1) m 4 



(82) 



meaning that, as could be expected, the contribution of 
the inflaton quantum fluctuations is significant only when 
inflation starts near the Planck scale, while it becomes 
negligible for smaller values of the initial energy density. 

We are now in a position where one can compute the 
new mean and variance exactly. As already mentioned, 
the new probability function is given by Eq. I|72f) with 
00 replaced by a. This immediately means that (Q) and 
0q are given by Eqs. if?^ and (j?4"|) with O replaced by 
0. The evolution of (Q) versus the number of e- folds 
is displayed in Fig. [SJ At initial times, one has (Q) ~ 
Q m and o 2 q ~ a 2 ~ 0g, the last approximate equality 
coming from the property established before Eq. (|82|l . 
On the other hand, the mean value and the variance of 
the quintessence field at late times read (if the initial 
condition for the inflaton field is large enough in order to 
reach the regime where a 3> Qi n ) 



(Q) 



(83) 



3n/2 



2n-2 



2n-2 



(81) 



I 

where in the above equations a is the a it) calculated at 
the end of inflation, which is given by 
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Wn 



3n(n + 2) m 4 



(n + 2)(7n + 4) V m 
3n(n + 1) to 4 



(84) 

where we have used Eq. I)82|l and an approximation for 00 
valid at late times. One notices that the two quantities 
of Eq. H83f) growing function of tp- m . This can be checked 
explicitly in Fig. where the final values of (Q) versus 
ifi n are represented. 

There is another important consequence that can be 
deduced from what has been discussed so far. As already 
mentioned, the energy density of the quintessence field at 
the beginning of the radiation era (i.e. at the end of in- 
flation) must be such that lO -113 ?™ 4 ^ < pq < 10~ 9 to 4 ; 
if one wants Q to be on tracks today. Assuming for sim- 
plicity that the field starts at rest, this implies 



Qmin _ ^q_107/c« < Q < 



Qn 



(85) 



Since, during inflation, Q is a stochastic quantity, one can 
study the probability for a given model (that is to say for 
a given choice of the power indices and of the initial con- 
ditions Qj n and </?i n ) that the field is in the appropriate 
range. This probability allows us to evaluate the likeli- 
hood of the various models under considerations, reject- 
ing those for which this quantity is small. In this way, we 
can thus constrain the value of the initial conditions for 
the fields and exclude a portion of the parameter space. 
Two remarks are in order here. Firstly, as is clear from 
Eq. (|72|l . the probability will depend on Q m . Therefore, 
although the late time evolution of (Q) is independent of 
Qin, this dependence is reintroduced via the calculation 
of the probabilities. Secondly, the result will also depend 
on (p- m . As a consequence, for a fixed value of n and a, 
one can hope to derive constraints on the initial value 
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of the quintessence field but also on the total number of 
e-folds during inflation. 

A given model will be accepted if a large part of its 
probability distribution calculated at the end of inflation 
is contained within the allowed range. A rough estimate 
of this constraint can be obtained by simply requiring 
that the mean value of the distribution falls in this range 
within one square root of the variance. Using the previ- 
ous results of this section, one finds that imposing that 
(Q) ± <tq falls between Qmin and <5max yields 




a < Q n 



(86) 



Therefore, in order to obtain a probability of order one 
must have that a < Q m ax- This implies an upper bound 



on (fi n that can be roughly estimated to (neglecting the 
contribution from inflaton fluctuations) 



< 



2/(r 



10 10(a-2)/[o(n+2)] _ /g^ 



It is easy to see from the above formula that this con- 
straint is quite stringent. Let us also notice that this 
also a constraint on the total number of e-folds during 
inflation, very roughly speaking N T < io 20 ( Q - 2 )/K"+ 2 )]. 

More precisely, if one uses the probability density func- 
tion given by Eq. (|72|l (with, as already discussed at 
length, do replaced by a if one wants to take into ac- 
count the inflaton noise at the end of inflation), one can 
calculate the exact probability for the quintessence to be 
on tracks today. One arrives at 



a, <j>hi,n) = 



1 



Erf 
Erf 



;max Q'm 

Qm'm H~~ Q'n 
V2<7 



Erf 



Qm'm Qh 

V2a 



Erf 



Qmax H~ Qi 

V5a 



The detailed behavior of this probability as a function 
of the initial conditions ip- in and Qin is shown in Fig. [7| 
for n — 2 and different power indices a. We see that 
a large portion of the parameter space can be excluded. 
Actually, for n = 2, the largest allowed initial condition 
is (fin ~ 7.6 x 10 5 m pl , corresponding to the Planck en- 
ergy density. However, at 99% of confidence level, initial 
values larger than 27m pi (a — 2), 1.3 x 10 3 to p; (a = 6) 
or 3.0 x 10 3 to p , (a — 11) are excluded. This also shows 
that, since the confidence region enlarges with the power 
index a, large values of a are statistically more favored 
than small values. 



V. DISCUSSION AND CONCLUSIONS 

In this section, we briefly discuss other aspects of the 
question studied in this article and present our conclu- 
sions. 

In the preceding sections, we have calculated the evolu- 
tion of the quintessence field during inflaton taking into 
account the quantum effects. What about these quan- 
tum effects in the subsequent cosmological eras? After 
the reheating, during the radiation and matter domi- 
nated phases, it is easy to show that the quintessence 
field evolves classically starting from the value reached at 
the end of inflation and determined by its random walk 
during inflation. Indeed, since the noise in the Langevin 
equation is controlled by the quantum fluctuations of the 
modes leaving the horizon, it rapidly becomes negligible 



as the accelerated expansion stops and the modes start 
to reenter the horizon. 

However, a new phase of accelerated expansion driven 
by the quintessence field is now taking place and, clearly, 
the above argument does not apply in this case. There- 
fore, one may wonder whether the influence of the quan- 
tum effects should not be taken into account when one 
computes the evolution of Q at present time. This could 
have important observational consequences, in particu- 
lar if the stochastic behavior of Q modifies the value of 
the equation of state. However, it is easy to demonstrate 
that this is not so. The problem is in fact very similar to 
calculating the evolution of the inflaton field during in- 
flation since, at present time, the quintessence field is no 
longer a test field but actually determines the evolution of 
the background. Therefore, even though the slow roll ap- 
proximation might not be so satisfactory in this case, we 
can at least estimate the relevance of these late stochas- 
tic fluctuations by simply setting n = —a into Eq. 11' I'D 
(since, in some sense, the Ratra-Peebles potential is noth- 
ing but an inflationary chaotic potential with a negative 
index). This equation, provided one substitutes tp c \ with 
Qd and V(tp c i) with W(Q c i), should be a reasonable es- 
timate of (<5Q 2 ) /Q 2 . Since the quintessence field is now 
on the attractor, its value must be of the order of the 
Planck mass. In addition, it is slowly rolling down its 
potential toward higher values. As a consequence, from 
Eq. H22|) . one gets the following rough estimate 



(SQ 2 ) 



Q 2 



< 



W(Q now ) 



10 



-123 



(89) 



18 



1 xicr 



8x10" 




Qin/mpi 



<l>in/ m PI 



2x10" 



- 


| >50% 




>io% : 


1 


> i% - 



20 30 

0in/ m PI 



40 



50 




i , 600 

tn/ m PI 800 



Qin/mpi 



0.10 



0.08 



0.06 



o 0.04 



0.02 



0.00 




1000 



500 1000 1500 2000 2500 

0in/ m PI 



a = 11 




1000 
4> in /m P , 2000 



Qin/mpi 




3000 



1000 2000 3000 4000 5000 6000 

0in/ m PI 



FIG. 7: Probability computed from Eq. I|88^ that the quintessence field is on tracks today for various models with n = 2 and 
different power indices a (from top to bottom, a — 2,6,11). For all the models under consideration, the maximum allowed 
value of the initial inflaton is (fii n = 7.6 x 10 5 m pi (coming from the fact that V\ n < Wit,). From the plots, one notices that values 
larger than 27m P! , 1.3 x 10 3 m pi and 3.0 x 10 3 m pi respectively (corresponding to an initial energy density Vi n of 1.3 x 10 _9 m pj , 
2.9 x 10" 6 m p , and 1.5 x l(T 5 m pi ) are excluded at 99% CL. One also notices that, since the confidence region enlarges for 
increasing values of a, large values of the power index a are statistically preferred. 
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meaning that any stochastic deviation from the classical 
trajectory is completely negligible at present time. 

Let us now end this work by reviewing what are the 
main conclusions of our study. The main result is that, 
taking into account the quantum effects during inflation 
is important since the stochastic diffusion term domi- 
nates over the classical drift term and that, as a conse- 
quence, the value of (Q) can differ from Q c \ by several 
orders of magnitude. For the first time to our knowledge, 
we have given an analytical estimate describing the evo- 
lution of (Q) during inflation, see Eq. 173fl . 

Another new result is the fact that, requiring the 
quintessence field to have a large probability to be on 
tracks today, allows us to put quite stringent constraints 
on the initial conditions Q m and ip- m . Typically, the 
quintessence field must start from small values. We have 
also established that large values of a are favored (we no- 
tice in passing that, if we consider the Ratra-Peebles po- 
tential only and not the SUGRA one, the same conclusion 
is reached from the constraints that exist on the equation 
of state today). Another interesting result is that the in- 
fiaton field must also start from quite small values. This 
implies that the total number of e-folds during inflation 
is also limited. On the other hand, we have remarked 
the existence of an attractor for (Q), see Eq. ifTfijl. due to 
the fact that the final value of (Q) is independent of Q ln . 
However, a dependence in the initial conditions is rein- 
troduced in the calculation of the probability which has 
allowed us to put the constraints mentioned just before. 

One of the main purpose of our paper was also to 
study the influence of the inflaton noise on the evolu- 
tion of the quintessence field. The approximation con- 
sisting in neglecting the inflaton fluctuations has been 
shown to be justified in most cases, basically because the 
corresponding contributions to the mean value and/or to 
the variance are proportional to Vin/m^ , see for instance 
Eq. (|83|l . Even in the extreme case of Planckian initial 
conditions for the inflaton field (i.e. V ln ~ frit .), the in- 
flaton noise is unlikely to modify (Q) /m rl by more than 
one order of magnitude compared to what is obtained 
taking into account the quintessence noise only. 

It is also interesting to compare these results to those 
obtained in the paper [l^ which was the first to take 
into account the quantum effects in the calculation of the 
evolution of (Q). Basically, our findings confirm and/or 
justify the results of Ref. [Hj and somewhat extend their 
validity. We have recovered the same equation for the 
variance and our new equation l|73|l for the mean value 
of Q confirms the conclusions that can be drawn from the 
figures of Ref. |15| , namely that the quantum effects can 
play an important role during inflation. In Ref. the 
inflaton noise has not been considered and, as mentioned 
above, we have demonstrated that this is, in most cases, 
a good approximation. 

Finally, let us describe some questions that are left 
unanswered and some possible improvements to the 
present study. In order to be able to find analytical so- 
lutions, we have modeled the classical drift term with a 



reflecting wall. The price to pay is that we have lost the 
dependence in the parameter a. The drift term acts dif- 
ferently for different Ratra-Peebles potentials while the 
wall repels the field regardless to a. Although we do not 
expect a strong dependence, it would be interesting to 
quantify this effect. The problem is that, if one includes 
the exact classical term, then the Langevin equation is no 
longer analytically solvable. The only way out seems to 
numerically integrate this equation. However, even this 
solution could be difficult because a term like Q~ a can 
rapidly become very large and, hence, problematic from 
the numerical point of view. 

Another interesting question would be to study what 
happens when one considers the case of a colored noise 
since it is clear that a white noise is not, physically, the 
most relevant case. Concretely, this amounts to replace 
the Heaviside function in the expansion of the field by 
a smooth function and, in principle, this could affect 
the evolution of the quintessence field during inflation 
although, again, we do not expect a very important ef- 
fect. 

For the moment, we postpone the study of all the issues 
to future works. 
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APPENDIX A: REFLECTING WALL FOR THE 
INFLATON 

Another way to look at the quantum evolution of the 
inflaton is the following. From Eq. i|31[l . it is clear that 
the variance of ^(t) is a growing function of time. How- 
ever, in order for the solution of Eq. (|28|) to be defined 
we need to impose the condition ^S>(t) < 1. Otherwise, 
this equation is clearly meaningless and, of course, if this 
condition is not satisfied, the series of Eq. is not 
convergent. This is probably the reason for the problems 
encountered before. The above condition can be thought 
of as constraining the random walk of \& with a reflect- 
ing wall. As explained in the preceding section, see also 
Ref. j2(| , the resulting probability distribution for ^ with 
the wall located at $ = 1 is given by 

-* 2 /(2(* 2 )) -(*-2) 2 /( 2 (* 2 )) 

P{% t) = 7=— , (Al) 

and the corresponding probability distribution for ip, ob- 
tained via the relation P[ip,t] = P[^(ip),t]\d^ /dip\, be- 
comes 



P{<p,t) = P[*(<p)] 



2 / (pc 



Vol \ V> 



(A2) 



This probability distribution has a finite mean value that 
can be expressed as [this equation can also be obtained 



20 



by using the link between \? and ip given by Eq. I|28(l and 
the probability distribution of given by Eq. QAljl ] 



e-V«* 2 )) , 

Vcl 7== — \/ 77-*- 1 



/4 



4(* 2 ) 



(A3) 



where I v {x) is the modified Bessel function of the first 
kind. For small values of ('J 2 ), the mean value reads 



(A4) 



yielding therefore the same results as (|33f) . Let also notice 
that all higher moments are divergent in accordance with 
the discussion presented in the section on the evolution 
of the inflaton field. 



APPENDIX B: PERTURB ATIVE SOLUTION 
FOR THE QUINTESSENCE FIELD 



Solution at First Order 



In this appendix, we present the explicit expressions 
(as a function of the classical inflaton field and/or of the 
number of e-folds) of (SQ 2 ) given by Eqs. (|6f J) and H62(l 
and of (6Q 2 ) given by Eq. We start with Eq. {o2jl 

which is the sum of two terms. In order to evaluate the 
first term coming from the inflaton noise, i.e. 
one must calculate the two-point correlation function of 
Sipi. Using the solution of Eq. I|19|l . one obtains 
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In the above expression, Q(z) is the Heaviside function, i.e. is zero if z < and one otherwise. Then, using the link 
between the cosmic time and the classical inflaton field in the slow-roll approximation, the remaining integrations can 
be easily performed since they just boil down to integrating power-law functions. One obtains 
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The evolution of (<5Q 2 ) |^ as a function of the number of e-folds is displayed in Fig. El The main feature of this 
formula is that it is proportional to (Wo/Vb) x (Wo/m*) which is, as discussed before, very small. This is confirmed 
by the plot, see the left panel in Fig. [SJ 

Let us now turn to the term sourced only by the quintessence noise. Using the fact that the quintessence noise is 
white the second term of Eq. I|tjl|) becomes 
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Then, using the expression of the inflaton and quintessence potentials, straightforward calculations lead to the following 
formula in terms of hypergeometric functions 
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FIG. 8: Left panel: evolution of (5Q\) /Q^i calculated in the perturbative regime for two different set of initial conditions (i.e. two 
different values of ipi n and Qi n ) and for a model with n — 2 and a = 6. For each set of initial conditions specified explicitly 
in the figure, the contribution coming from the quintessence noise and the contribution originating from the inflaton noise are 
represented. The perturbative treatment is under control only in the case tpi n — 221m p , and Q m — 0.06m pi corresponding to 
initial energy densities given by Vi n = 3.7 x 10 _8 m pi and Win = 2.8 x 10 _114 m p , (since all the corresponding curves lie in the 
region (<5Qi)/Qci < !)■ The plot demonstrates clearly that, in the perturbative regime, the inflaton noise is totally negligible. 
Right panel: evolution of {5Q2)/Q C \ for the initial conditions <pi n = 221m pi and Qi n = 0.06m P! . The five curves correspond 
to the five contributions of Eq. 1641 . Again, it is obvious that the four contributions originating from the inflaton noise are 
unimportant. 



where we have used the definitions 
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and where we have assumed that n ^ 2. The number /i is always positive for n < 6 which is the case we are mostly 
interested in and the number v is, on the contrary, always negative. The evolution of {SQf) is represented in Fig. EI 
The presence of the hypergeometric functions is linked to the fact that the classical quintessence held is not totally 
frozen. Looking at Eq. I|5U|) . one sees that the term responsible for the slight evolution of Q c \ is in fact proportional 
to A. Indeed, in the limit A — > 0, or equivalently T — > 0, the classical quintessence held becomes Q = Q ln , i.e. is 
actually frozen. In this case, one has 
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and Eq. (|B4f) reads 
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which is exactly the expression derived in Ref. 0]. This term is proportional to the factor Vo/m pi which is much 
larger than (Wq/Vq) x (Wo/mt,)- As a result, we expect the term originating from the quintessence noise to dominate 
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over the term originating from the inflaton noise and, as already mentioned, this conclusion is confirmed by the plots 
in Fig. EI A priori this conclusion is valid only in the regime where the above expressions have been established, i.e. in 
the perturbative regime. 

We have noticed before that the solution expressed in terms of the hypergeometric function is valid provided n =/= 2. 
The case n = 2, for which the evolution of the classical quintessence field is given by Eq. 1)5 requires a special 
treatment. In this case, one finds 
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where "f((3,x) is the incomplete gamma function defined by 

709, a:) = / e'H^^dt. 
In the limit T — > 0, using the formula 7(/3, a;) ~ 
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which is again the expression found in Ref. |15|. specialized to the case n = 2. The previous conclusion, namely that 
this term dominates the contribution coming from the inflaton noise, is not modify in this particular case. 



2. Solution at Second Order 



We now turn to the calculation of the mean value of 5Q2 given by the sum of five contributions as can be seen in 
Eq. (|64|l . In principle, the calculations of these five terms can be performed in the general case, where the classical 
quintessence field in given either by Eq. (|50|l or Eq. I|51l) . However, for simplicity, we give only the expressions 
corresponding to the case where Q c \ is frozen, i.e. to the limit T — ► 0. In this case, all the integrations become trivial 



since they only involve integrals of power-law functions. The first term, {5Q%) 
It is given by 
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originates from the term 5<p2 only. 
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The main feature of the above term is the presence of the overall factor Wo/m^ x (<5i n /™ P! ) a 1 (of course, the 
powers of the inflaton field can also play a role but for a crude order of magnitude estimate, one can ignore them) . 
The second term participating to the expression of (6Q2) comes from 8ip\. It is given by 
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We see that this term also scales as Wo/iri^ x (Qi n /^ P! ) a ■ Therefore, we expect the previous contribution and 
this term to be of the same order of magnitude. This can be checked explicitly in Fig. [S] The third term comes from 
dfiSQi where it should be understood that, in SQi, only the term proportional to the inflaton noise matters. One 
obtains 
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In comparison with the two previous terms, we see that there is an extra overall factor equal to Wq /Vq x (Qi n /?7i pi ) 
As a consequence, this term is expected to be sub-dominant with respect to the two previous contributions and, again, 
this can be checked explicitly in Fig. El Of course, one could try to compensate the smallness of Wq/Vq by a large 



value of (Qi n /m Fl ) a 2 , i.e. by a large value of the index a but this would lead to very artificial, hence unphysical, 
models. Then, the fourth term, originating from the term 5Q\{£,q), can be written as 
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If we compare this term with the two first contributions, we see that there is an extra overall factor equal to 
(<2i n /m pi )~ 2 . Since Qi n /m< F i is small one expects the previous contribution to be dominant. In Fig. [HI this con- 
tribution is represented by the solid black line which is indeed the most important one. Finally, the last term coming, 
from <5Q 2 (£), reads 



(6Q2 



2a 3 {a + l){a + 2) W W 2 f Q ia 

i-n/2 



3n 3 
32 



i 4 Vr? 

PI u 



n(4 — 3n)(8 — n) V m p 
4 



1 



2n(n- 1) 



Hi 



n[n 



2)(4- 



<£it 
rn 



'fid 

i—n 



-3a-5 



2-2n 



2(3-n)(4-n)(8-n) 



m D1 



2-3ri/2 



fin 

2-n 



2-2n 



(fir. 



2-n 



2-3ri/2 



(B15) 



This term is proportional to {Wq/Vq) 2 which is a tiny 
factor. Therefore, the above contribution is expected to 



be the smallest contribution to {8Q2) and, in fact, to be 
totally negligible. This is confirmed in Fig. [SJ 
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